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THE SOLUTIONS OF NON-HOMOGENEOUS LINEAR DIFFERENCE 
EQUATIONS AND THEIR ASYMPTOTIC FORM* 

BY 

K. P. WILLIAMS 

Introduction. 

Very important progress has recently been made in the analytic theory of 
homogeneous linear difference equations. 

Galbeun f has used the Laplace transformation to derive important ex- 
istence theorems, and has investigated the nature of certain principal solu- 
tions for large values of the variable. 

At about the same time Noelund J applied the theory of factorial series 
to linear difference equations, showed the existence of the same solutions, 
and gave their asymptotic form. In addition to the case of polynomial co- 
efficients he has considered the case where the coefficients can be expressed 
in factorial series. In a more recent paper he has investigated an equation 
with still more general coefficients. § 

By means of a method of successive approximation, and a suitable extension 
of a contour integral due to Guichard, Carmichael || has independently 
shown the existence of solutions, and has found the bounds of increase and 
decrease of the solutions in a direction parallel to the real axis. 

The latest important contribution to the general theory is by BiRKHOFF.If 
He employs a matrix notation, and shows in a direct manner the existence 
of certain intermediate solutions and of the principal solutions. The asymptotic 
form of these solutions is determined throughout the complex plane. A 
modification of the integral used by Carmichael plays an important role. 

* Presented to the Society February 24, 1912. 

jComptes Rendus, vol. 148 (1909), p. 905; vol. 149 (1909), p. 1046; vol. 150 (1910), 
p. 206; vol. 151 (1910), p. 1114; Acta Mathematica, vol. 36 (1912), pp. 1-68. 

tComptos Rendus, vol. 149 (1909), p. 841 ; Bidrag til de lineaere IHffemsligningers 
Theori. Dissertation, Copenhagen, 1910. 

§M6moires de I'Acad^mie Royale des Sciences et des Lettres 
de Dancmark, Copenhagen; 7th series, vol. 6 (1911), pp. 309-327. 

|| These Transactions, vol. 12 (1911), pp. 99-134. 

"[These Transactions, vol. 12 (1911), pp. 243-284. 
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He finds also the form of the periodic functions which connect the principal 
solutions, and formulates for the first time the properties which are sufficient 
to characterize a given system. 

While the analytic theory of homogeneous linear difference equations has 
thus been extensively treated, no general theory has been developed for non- 
homogeneous equations,* although a number of equations of particular form 
have been considered (see Carmichael, loc. cit., Introduction). 

The purpose of the present paper is to prove the existence of certain principal 
solutions of the system of non-homogeneous equations 

n 

gi(x+l) = ^2a ij (x)g } (x) (i = l, 2, •••, n -1), 

(1) 

n 

g n (x+ 1) = ?L,a n j(x)g j (x) + b(x+ 1) , 

3=1 

and to develop their asymptotic form. The functions a t j (x) and b (x) will 
be assumed to be rational,! with poles of order at most ju and v respectively 
at infinity, so that we can Write 

a,ij{x) = x» f ay + ~ + • • • J , 



.(*) = *"(& + ^ +..-), 



for I x I > B . 



On account of the gain in simplicity of the formulae we have taken a system 
where all but one of the equations are homogeneous. It is, furthermore, evident 
that the system (1) would be no more general if all the equations were non- 
homogeneous, since the solution of such a system could be obtained by adding 
together the solutions of the n systems, of the above type, which would be 
obtained by striking out the non-homogeneous term from all the equations 
except the first, second, . . . , last, in turn. Thus the theorems for (1) can 
be extended readily to a system of n equations, all non-homogeneous. 

The method employed is based on a 'method analogous to that of varia- 
tion of constants in the theory of linear differential equations. This method 
has previously been supposed to yield only formal results.J It will appear, 

* It is possible to reduce a non-homogeneous equation to a homogeneous equation. The 
application of the general results for a homogeneous equation will show the existence of solutions, 
but gives no direct means of studying their properties. This method is suggested by Wal- 
lenberg and Guldberg, Theorie der linearen Dijferenzengleichungen, p. 87. 

t Most of the results which will be obtained apply with slight modification to the case where 
an ( x ) and b ( x ) have the form of rational functions only at infinity. 

t Wallenberg und Guldberg, p. vii, "... die Methode der Variation der Konstanten, 
die wegen der mit ihr verbundenen Summationen meist nur formale Bedeutung hat . . . ." 
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however, that in general, if not in all cases, the sum formulae that arise can be 
interpreted so as to give solutions analytic in the finite plane except for 
poles. 

There is a difference of treatment according as ju > , ju < , or 
fi= 0. 

In § 2 the reduction of the problem to that of making valid one of four 
formal types of sums is effected. 

In § 3, under the hypothesis that n > , a direct summation to the right 
is used to derive a first principal solution of (1), which we denote by <7n ( x ) , 
</2i ( x ) , • • • , g n \ ( x ) . This solution is analytic in the finite plane save for 
poles at determinate points near the negative real axis. In § 4 it is proved 
that in the right half plane this solution is asymptotic to the series formally 
satisfying (1). 

It is shown in § 5 that the formal solution of (1) obtained by a direct sum- 
mation to the right, such as is used in § 3, is independent of the fundamental 
set of solutions of the homogeneous system associated with (1) which is used 
to determine it. This principle is used in § 6 to determine the asymptotic 
form of the first solution in the left half plane. 

A summation by means of a contour integral to the left is employed in § 7 
to determine a second principal solution, gu (x) , gm (x) , •••, g n2 (x) , 
analytic save for poles at determinate points near the positive real axis. This 
solution is asymptotic to the formal series in the left half plane. 

In § 8 the relation between the principal solutions is determined. 

It is seen in § 9 that, if n < , the same theory holds except that the 
roles of the right and left half planes are interchanged. Similar results are 
proved in § 10 for the intermediate case jw. = 0, except in certain extreme 
cases. 

An application of the general theory is made in § 11 to a single equation 
of the nth order. 

The writer wishes to express his appreciation and thanks to Professor 
Birkhoff for the valuable suggestions and aid that have been received from 
him. 

§ 1. The fundamental existence theorems for a homogeneous system. 
The system 
(2) 9i(x+l) = T,a ij (x)g 3 -(x) (i = 1, 2, ••• , n) 

is the homogeneous system associated with (1). As we shall make use of its 
solutions, it will be convenient to enumerate some of their properties. We 
state below the existence theorems for (2) as they are given by Birkhoff. 
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There exist in general n sets of series: 

a u (x) = x»* ( Pi e-n x x ri {s u + *-% +•••)> 

(.(i) \ 



«„(*) = a:- (p y *-*)-a: r >^ + -| +•••], 

each of which will formally satisfy (2), but which in general diverge. 
The quantities pj are the roots of the determinant characteristic equation 

(4) lay- 8,7 p | = 0, 

where 5y = for i 4= j , and 5,-,- = 1 , and they are so ordered that 

(5) | Pi | 5 | Pi | !> • • ■ % | P» | • 
Furthermore, they are all supposed different from zero, that is 

(6) A = | aij | = pi pi • • ■ p n 4= . 
It is assumed that such series exist, and also that 

(7) s=KI*o. 

We can now state the theorems we shall use. 

There are two fundamental sets of solutions of (2), the first and second principal 
solutions;* they may be denoted by hu (x) , h^ (x) , •••, h„j (x) (j = 1, 
2, ••• , n), and hi, (x) , hy (x) , ••• , h n} - (x) (j = 1, 2, ••• , n) , re- 
spectively. They hate the defining property that in any right half plane 
h» (x) ~ Sij (x) ,f while in any left half plane ha (x) ~ sy (x) (i, j = 1 , 
2 , • ■ • , n ) . TAe functions ha (x) are analytic in the finite plane save for poles 
at the poles of oy ( x ) , and tf Ae zeros of the determinant | ay ( x ) | , or points 
congruent% to these points on the left; while the functions Ay (x) are analytic 
save for poles at the poles of ay (x) , and points congruent on the right. 

While all the results which will be derived could be obtained by making 
use of only the principal solutions of (2), it will simplify the work to make 

* Our first and second principal solutions are the second and first principal solutions, respec- 
tively, as given by Biekhopp. For the present purpose it is convenient to take them as above. 

t This means that for any c the difference between ha (x) and A; + 1 terms of s,-,- ( x ) can 
be written x* x (p,e-" )*£'>• ft (a; )/:!*, where lim f* (x) =0, uniformly, for — Jir^arg (x — c) 

Z=QO 

^ §t or, in a special case, for — $?r — t ^ arg (x — c)^Jir — <= , e arbitrarily small (see 
Galbrun, loc. cit.). We say then that hn(x) is asymptotic to s« (x) , or in symbols 
ha (x) ~Sij (x) , and call f * ( x ) the remainder coefficient. 

t The points ••• x — 2 ,x — \ ,x ,x + \ ,x +2, ••• axe called a congruent set of points. 
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use of the properties of one of the two sets of intermediate solutions that are 
given by Birkhoff (loc. cit. § 3). 

There exists a set of intermediate solutions h[j (x) , h' 2 j (x) , • • • , h' n j (x) 
(j = 1,2, • • • , n) , defined only at a sufficient distance from the axis of reals. 
The elements of these solutions are analytic in the region where they are defined, 
and are such that h'ij (x) ~ sy (a;) with regard to x in the left half plane* and 
with regard to v (x = « + W ' — 1 ) , in any right half plane. f 

Although these intermediate solutions are not defined near the real axis, 
they will be of service to us on account of the fact that they maintain their 
asymptotic form with regard to v in any right half plane, and because the 
periodic functions which connect them and the second principal solutions are 
of a simple nature. The relation between these intermediate and the second 
principal solutions is given by the formula 

hj (x) = <pij (x) h'n (x) + <pij{x) h' i2 (x) + ■•■ 

(8) 

+ <pj-i,j (a;) h, j-i(x) + ha (x) (i, j = 1, 2, ••• , n) , 

where <pu ( x ) , <p 2 j ( % ) , • • • , <Pj-i, j ( x ) are periodic of period 1 , and are 
analytic at a sufficient distance from the real axis.f 

When we put j = 1 in (8) we see that the solution hn(x) , • • • , h nl ( x ) 
(i. e., the one corresponding to the largest \ p\) maintains its asymptotic 
form in the right half plane with regard to v . These functions form the ele- 
ments of the first column of the determinant \hij(x)\. It is also proved 
in Birkhoff 's paper that | ha ( x ) [ , or any minor formed from its first k 
columns ( k < n ) , is asymptotic to |*# ( x ) | , or its corresponding minor, 
with regard to » in the right half plane. 

There exists a second set of intermediate solutions h'ij ( x ) , such that 
h'n (x) ~ Sij (x) with regard to x in the right and with regard to v in the left 
half plane. 

§ 2. The formal solution and the operation of summation. 

We are thus assured of the existence of solutions of a homogeneous system. 
This will enable us to determine easily a formal solution of (1). For this 
purpose let Ijj (x) , hj (x) , • ■ • , l n j (x) (j = 1 , 2, • • • , n) be any funda- 
mental system of solutions of (2). Write 

n 

(9) gi(x) = ^,o>j (x)kj (x) (i = 1, 2, ••• , n), 



* That is, asymptotic in the sense already defined. 

t This means that the difference between hij ( x ) and k + 1 terms of s,-,- ( x ) can be written 
x» x (pje-*) x x'jZi(x)v- k where lim fi(x) = 0, —\v= arg (x — c)^£x. 

r=oo 

% Birkhoff, loc. cit., p. 266. 

Trans. Am. Math. Soc. 14 
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and let us determine the quantities wi|(a;) , o>2 (x) , • • • , w,; (a;) so that gi(x) , 
92 (x) , • • • , g n ( x ) will be a solution of the non-homogeneous system. If we 
substitute the latter expressions in (1), after first putting 

o>j ( x + 1 ) = wj ( x ) + Awy ( x ) , 

we find in virtue of (2) 

n 

J^Aoij (x) kj{x+ 1) = (i = 1, 2 ,[-.., n-1), 

ZA Wy (a;)^(a!+l) = 6(a;+ D- 

This" is a system of n linear non-homogeneous equations in the n quantities 
Afa>i ( x ) , A&>2 ( a; ) , • • • , Aw„ ( a; ) . Let i ( x ) represent the determinant 
\kj (x) | . Further, let iy ( a; ) be the minor of the elenient in the last row 
and jth column of L ( x ) . From (10) we have 

(11) A. J (x)=(-l)^ bix \\l L i_[ X ) +1 \ 

We are thus naturally led to the consideration of the simple non-homogeneous 
equation 

(12) Af(x)mf(z+l)-f(x) = 0(x)* 

There are two series which formally satisfy (12), namely, 

(13) f(x) = -6(x)-0(x+l) , 

(13') f(x) = d(x-l) + 0(x-2)+ ■■■, 

which, if they converge, will be actual solutions. These will be called the 
direct sums to the right and left respectively. 

Besides the direct sums, a formal solution of (12) exists in the form of a 
contour integral due to Guichard.f We shall use it in the form adopted by 
Carmichael.J Let 6 (x) be analytic outside of the region D , extending to the 
right and lying within a finite distance of the positive half of the real axis 
(see fig. 1); then if the contour L = oo AB oo is drawn so as to enclose D, 
and passes between x— 1 and a; in the manner shown in the figure, the function 

* The solution written in the form ^6 (x) is called a sum formula where the operation S is 
such that S«(i)-J«(i)=((i!), 

X+l X 

t Wallenberg und Guldberg, loc. cit., p. 12. 
t Loc. cit., p. 119. 
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formally satisfies (12) by Cauchy's formula, since / ( x + 1) — f (x) reduces 
to the integral around ABB' A' , and the only singularity of the integrand 
within this region is at t = x* 

If now 8 (t) decreases along A oo and B °o in such a way that the integral 
converges, (14) will furnish an actual solution of (12). 

It is obvious that in case the region D extends to the left, we can take a 
contour L extending to the left so as to again enclose D . The sign of the 
integral must in this case be changed if we integrate so as to keep the area 
within L on our left. 



d-i 




Fig. 1. 



We therefore have four formal solutions of an equation of the type (12), 
the direct sums to the right and left respectively, and contour integrals about 
contours extending to the right or left. 

When one of the above sums is applied to (11) we see that (9) will furnish 
a formal solution of (1). Our main ■problem therefore is to determine whether we 
can so choose the system of solutions h, ( x ) , ly ( x ) , • • ■ , l nJ - ( x ) ( j = 1,2, 

• • • , n) ,of (2), that the application of one or more of the four sums to (11) will 
lead to an analytic solution of that equation and therefore of (1). 

§ 3. The first principal solution for n > . 

We choose first as the solutions of (2) the first principal solutions, hjj (x) , 
h-u ( x ) , • • • , hnj ( x ) (j = 1 , 2 , • • • , n ) . Take, further, as the solution of 
(11) the direct sum to the right. Thus we have formally 

(15) t o-(a ; ) = (-l)^[ b(a; + 1) ^ (a: + 1) + 6(a:+2) ^ (a:+2) +---1 

where H ( x ) is the determinant | hjj ( x ) | , and Hj (x) is the minor of this 
with reference to the element in the last row and jth column. 

* In case 6 ( x ) is analytic within D , the integral in (14) reduces to (13) by a familiar theorem 
in the theory of residues, since the only singularities of the integrand are &tx,x + l,x + 2, 

• • • . That is, we are led to (14) on writing (13) as a series of residues. 
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Each of the terms of (15) will be analytic save at the zeros of H(x) , and 
the poles of b ( x + 1 ) and H } - ( x + 1 ) , or points congruent to these points on 
the left. 

As has been stated, hij ( x ) is asymptotic to s,j ( x ) in the right half plane, 
so that H ( x ) is asymptotically represented by the determinant \s { j ( x ) | . 
From (7) it follows that H ( x ) will have no zeros in the right half plane at a 
sufficient distance from the origin. Since the quantities ha (x) form a 
system of solutions of the homogeneous system, we see easily that H ( x + 1 ) 
is equal toA(x)H(x), where A ( x ) is the determinant \a i} - ( x ) | . Therefore 
we have H (x — 1) = H (x) / A (x — 1), and since ^4 ( a; ) has no singu- 
larities in the finite part of the plane exterior to a circle of radius R , it follows 
that H ( x ) will not vanish in the finite plane at a sufficient distance from the 
negative real axis, say at a distance greater than R'. 

From the situation of the singularities of hy (x) it is apparent that Hj (x) 
is analytic in the finite plane at a distance greater than some quantity R from 
the negative real axis, while b ( x + 1 ) is analytic if x -f- 1 lies outside a circle 
of radius R . 

Let T be the greatest of the quantities R, R' , and R. About the origin 
draw a circle of radius T , and the half -lines tangent to it and parallel to the neg- 




Fig. 2. 

ative real axis (fig. 2). Denote by P the region enclosed by these lines and the 
semi-circumference in the right half plane. Then if x is any point in the finite 
plane exterior to P , all the terms of (15) are analytic. If x is within the region 
P , some of the first terms may have singularities, but all the terms from the 
mth remain analytic, if m be so chosen that x + m is exterior to P . 

It remains now to discuss the convergence of (15). Consider the mth term. 
We will show that it is less in absolute value than the mth term of a rapidly 
converging series. To do this, replace the elements of Hj ( x + m ) and H ( x+m) 
by their asymptotic forms. When this is done it is evident that we can remove 
the factor (a: + m y {x+m)+r > ( Pk e ~*) x+m from the Arth column (k = 1, 2, 
■ ■ • , n) of H (x -\- m) . All these factors can be removed from Hj(x + m) 
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except the one for k = j , so that we have* 



b(x + m)Hj(x + m) (x + m)" I df \ 

H(x+m) ~ (x+m)» ( -* im '> +r '( Pj e-») 1 * m \ ^ x+m^ /' 



(16) 

where 

(16') d, = b -f : 



Sj being the minor of the element in the last row and jth column of S = \ s i} - \ . 
If we put x + m = mi , we see that the typical term of our series can be 
written 

m[C mi 



m 



^"•i 



Minn), 



where c and C are constants, and M ( m x ) is a uniformly limited function 
for m large and x in the vicinity of any point Xi . Since mi increases indefi- 
nitely in the successive terms of the series, it is evident that the above 
expression is the general term of a uniformly and rapidly converging series, 
irrespective of the values of c and C , provided only that n > , which we are 
assuming for the present. 

By Weierstrass's theorem the series (15) will therefore represent an analytic 
function, provided x is exterior to P; for this makes all the terms analytic. 
If x is within P, the first terms of (15) may, as seen above, have poles, the 
remaining part of the series converging uniformly, and accordingly «,• (x) 
may have a pole. 

The values of wi ( x ) , • • • , w» ( x ) when used in (9) will furnish a solution 
of (1) which is analytic in the finite part of the plane save for poles within P . 

The solution which is obtained in this manner can be extended to the left 
as follows: let Aj% (a;) be the minor of the element in the jth row and ith 
column of A ( x ) , and solve ( 1 ) f or g\ ( x ) , • • • , g n ( x ) , as we may do since 
A ( x ) ^ . This changes the equations to the form 

d7) 9i (x) = ±(-iy«^0t g . (x+1) + { _ 1)n+i+ i40± Hx+1) 

(i = l, 2, •••, n). 

The solution can now by (17) be analytically continued to the left. 

From this form of (1) we can also determine the location of the singularities 
of the solution without reference to those of hi,- ( x ) . Since the g's are analytic 
in the finite plane exterior to P , it follows that, in their continuation to the 
left, singularities will be introduced only within P at the singularities of 
b ( x + 1 ) , at the minors of A ( x ) , and at the zeros of A ( x ) . A singularity 
will in general recur at points congruent on the left. 

* The series asymptotically representing the sum, product, or quotient of any number 
of functions is found by adding, multiplying or dividing their respective asymptotic series 
according to the ordinary rules for convergent series, Borel, Les S tries Divergenles, p. 30. 
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By using the first principal solutions of the homogeneous system, and the direct 
sum to the right, we have thus shown the existence of a first principal solution 
0H (x) > 921 (x) , • • • , g n i (x) ,of (1), which is analytic in the finite plane, save 
for poles near the negative real axis, congruent on the left to the singularities of 
aij ( x ) , b ( x + 1 ), and the zeros of A ( x ) . 

§ 4. The asymptotic form of the solution in the right half plane. 

We shall now investigate the asymptotic form of the solution <7u ( x ) , 
02i ( x ) , • •• , <7„i ( x ) . Manifestly for this purpose it is merely necessary 
to consider more closely the terms of (15). The asymptotic form of the typical 
term has already been determined, and is given in (16). 

Before proceeding it will be convenient to introduce a simplification in 
notation. When / ( x ) is asymptotically represented by a power series s ( x ) 
of negative integral powers of x, instead of replacing / ( x ) by the expression 
which is obtained on breaking off the series s ( x ) at some point and inserting 
the remainder coefficient, (see footnote, § 1), we shall put 

f (x) = [a ; x], 

giving in each case the constant term a and the argument x in the series.* 
The leading factor of (16) can be written 



(18) — - - 



x liX (p j e-»)*x r > +li , n / , mV + ' l(l+m) 



1 + -) (p;e~*)" 



the first factor of which is independent of m . We shall omit, for the present, 
this factor, so that we have to consider the series 

» ( e M / p • ) m 
( 19 ) A ? J xrj -H-M(x+m) K-5 x+m]. 



-' ('+;)' 



Consider in the first place the first k terms of (19), k < \x\, that is, the 
finite series 

\ l + x) *V + ; 

(20) 

, | (e»/ P j) k 

~r ■' " ' ~r / t,\ r,— v+ „.( x+k ) [aj, x-f- k j, 

^(l + l) 

* This notation is similar to that used by Birkhopp, these Transactions, vol. 9 
(1908), pp. 373-395. 
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and let us examine the first term in detail. We have 

e* / Pi _ e*-{'>—+* <*H) } io« u+i/s) _ e ^-^ + - 1 Xi X_ 2 

~r _ ~r _« • • • • 



or 



i-H-dfc+U py py Pj X ' X 2 



If then, in the series dj + dj'V (a; + 1 ) + • • • , we put 

z+1 a:\ a;^ /' 
the first term of (20) can be written 

p y ~ t ~a;~ , "" i ~ ■«"* 

where E^ (a;) is a uniformly bounded function for \x\ sufficiently large and x 
in the right half plane, since it is composed of remainder terms in uniformly 
convergent series and the remainder coefficients in asymptotic series.* The 
terms of (20) from the second on can be treated in the same way; but we notice 
that their expansions begin with ar" , af~ 2 " , • • • , x~ iir ~ V)>l , respectively. The 
first k terms can therefore be written 

dj , h , l± , , huk+ Hah ( x ^ 

Pi + x + x 2 ^ h x» k ' ' 

where rj^k (x) is limited for x large in the right half plane. 

It remains now to consider the portion of (19) from the term m = k + 1 
on. We will find an upper limit for the sum of these terms. In the first 
place we notice that the factor 

' + ?) 

which occurs in the denominator, is greater in absolute value than some quan- 
tity K, and in general increases with m, since x is in the right half plane. 
Also the expression [df, x + m] can be* replaced by a quantity larger in ab- 
solute value. Consequently the portion of our series from the (k + l)th 
term on is less in absolute value than the series whose general term is 



M 



K I x |*<*-« 



* To prove that / ( x ) ~ a + ai / x + a 2 lx 2 + • • • , we need show only that ft ( x ) in 
f (x) = a,o + ai/ x + ••• + [a,h + ?k(x)] I x k remains bounded for every k , since 
ft ( x ) = [ aic+i + f t+i ( x ) ] / x , and hence lim ft ( x ) = . This principle will lje used in 

3 = 00 

what follows. 
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where M does not increase with \x\ and K does not decrease. This is a geo- 
metrical progression, whose ratio can be made as small as we please by taking 
| a; | large. The sum of the comparison series is therefore of the order of its 
first term, which is 1/ | x |"*, since we are starting with m equal to h + 1 . 
This is of the same order as the last term of the expression found above for 
the first m terms, and consequently (19) can be written [djj pj-, x], where the 
series can be taken to as many terms as desired, since k was arbitrarily large. 

An expression for ccj (x) will be obtained by multiplying [djj p,; x] by 
( — 1 ) n +J' +1 , (see (15)), and by the first factor in (18), which was removed from 
(19). 

The quantity we wish to determine is wy ( x ) h^ ( x ) . Therefore, since 
ha ( * ) = *"*( Pi e ~* Y x Tj [ s a> x], we have 

(21) ioj(x) hij (x) = x"-» [ ( - 1 )"+' +1 d, Sij I Pj ;x]. 
Letting j = 1,2, • • • , n and adding one finds 

(22) gn{x) = x v - lx [h i ;x] (i - 1, 2, •••, n) . 
where 

Jkl - ( _i ) .(*L!»_*i« + ... + ( _i ) ^£&y 
\ pi 9^ p» / 

From (6) and (16') we obtain 

(23) hi = — gj — ? ^ — - 1 ) i+1 PlP2 '" M- 1 Pi+1 '" Pn s v & 

(i = 1, 2, ••• , n). 

Let us now try to find series which will formally satisfy (1). Write 

(24) 9i(x) = x v -»(g i + 9 -^+ •••) (*-l,2, ■-., n>, 

and substitute in (1). When we develop the resulting equations according to 
decreasing powers of x it is evident, on remembering that ft > , that we can 
uniquely determine the successive constants g it g^ , • • • by comparing coef- 
ficients. The developments in all the equations start with x"; equating its 
coefficients we find 

n 

= 2 a H gj (i = 1 , 2 , • • •, n - 1 ) , 

n 

= Ha^+b. 

From this we have for the first coefficient in gi ( x ) 

hA ■ 

(25) gi= (-D n+i+1 -r- «-l,2, ■■-, n), 
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where A n % is the minor of the determinant A = | a»y | with reference to the ele- 
ment in the last row and ith column. 

We call the series g\ (x) , • • • , g n (x) the formal solution of (1). These 
series are obtained by direct substitution, and in general diverge. 

Since the expansions in the expression for gn (x) and in the formal series 
g, (x) (i = 1,2, • • • ,n) are of the same form, the question naturally arises 
as to whether their coefficients are the same. We can show that this is the 
case in the following way. We know that gn{x) , • • • , g ni ( x ) as given in (22) 
satisfy (1), and if we substitute this form of the solutions in (1) and compare 
coefficients of quantities of the same order we see that we must have precisely 
the same equations satisfied by h t , h^ , • • • that were satisfied by g it g { P , • • • 
( i = 1,2, • • • , n) . Since these equations had unique solutions, this means 
that corresponding coefficients in the series ga(x) and g, (x) (i = 1,2, • • • ,n) 
are equal. 

Thus we see that the first principal solution is asymptotically represented 
by the formal solution in the right half plane. 

Although we have shown easily that the constants in the series for gn (x) 
must be the same as those in the formal series </,- ( x ) , it is interesting to note 
that we can show by direct computation that hi is identical with gt . To do 
this we must eliminate the quantities py and s»y from (23). 

When we substitute the formal series (3) in (2) and compare coefficients of 
x* +ri we find 

n 

2(a»fc — 8ikP 3 ) s h j =0 (i = l, 2, ■■• , n), 

where 8«. = if i + k , and 5»» = 1 . Since py is a solution of "(4), the last equa- 
tions have solutions other than Skj = . The first n — 1 of them give 

n 

PjSij = ^,aikSkj (i = 1, 2, ••• ,n - 1)J. 

If now we introduce the quantities pi , • • • , py_i , py+i , • • • , p n into the succes- 
sive columns of the determinant &»y and eliminate the p's by means of the 
above expressions, we obtain the relation 

(26) Pl • • • Py_i Py+1 • • • Pn Snj = A n i Slj + A ni Sij + • • • + Ann Snj, 

where An and £»y are the minors of A and S respectively with reference to the 
element in the ith row and jth column, and where we have put Snj = Sj, 
the two being identical by definition. 

If in (26) we let ,; = 1,2, • • • , n , and substitute in (23), we find for i = 1 

(-1)"6 ^ i (-D"&Ai 
«i= — gj — 2-,A nh 2-.\— iy +1 si } Skj= -j = g\, 
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since 



U { L) Sl} * k > 10 if A + 1. 



In a similar way we can show h, = g t (i = 2 , 3 , • • • , n ) . 

We summarize the results of this and the preceding section in the following: 

Theoeem I. If n > , the system of non-homogeneous linear difference equa- 
tions (1) has a first principal solution gu ( x ) , g 2 i ( x ) , • • • , </„i ( x ) . All the 
constituent elements of this solution are analytic throughout the finite plane, save for 
poles at the poles of the functions a,-,- (x) ,b (x + 1 ) , at the zeros of the determi- 
nant \aij(x)\, and at points congruent to these points on the left. There exist 
n series gi(x) , g 2 (x) , • • • , g n (x) which formally satisfy the system (1), and 
whose coefficients can be determined by direct substitution. These series in gen- 
eral diverge, but when x approaches infinity in any right half plane, gn (x) is 
represented asymptotically by gi (x) . 

This theorem is supplemented by Theorem II of § 6. 

§ 5. A formal property of the direct sum. 

The solution of (1) which we have found was determined by making use of 
one of the sets of principal solutions of (2). The question now arises as to 
whether a different solution of (1) would be obtained if we used a different 
fundamental system of solutions of (2) and employed again the direct sum 
to the right. The answer is in the negative as will be shown. 

Let lif ( x ) , hj ( x ) , • • • , l„j ( x ) ( j = 1 , 2 , • • • , n) be any other funda- 
mental system of solutions of (2). From the theory of homogeneous systems 
we have 

n 

hi(x) = 2 <Pa (x) hi (x) , 
(27) hi (x) = J2 <Pa (x) h 2i (x) , 0' = 1, 2, ••-, n), 



Inj (x) = 2 <Pa ( x ) hi (x), 

i=l 

where the functions <pa ( x ) are periodic of period 1 , and the determinant 
$(x) = |*«,(a;)|#0. 

Denote by </,• ( x ) ( i = 1 , 2 , ■ ■ • , n) the solution of (1) arising from (27) . 
Then one obtains the formula analogous to (9), 

n 

(28) ?.-(«) = I>y (a;) l i} - (x) (i - 1, 2, •••, n), 

1=1 
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where the functions dj ( x ) ( j = 1 , 2 , • • • , n ) are to be determined. As 
before, let L ( x ) denote the determinant \ly (x) | , and L, ( x ) the minor of 
this with reference to the element in the last row and j'th column; then we 
have by (15), since the work to that point was purely formal, 



(29) 



Mx)-(-.^[ t( ' + t ^ + ^ + -] «-.... •-..■ 



From (27) it follows that Z ( x ) = H (x) •*(»). Similarly, we can show 
that 

4- (a) = *y (a) Ht (x) + * 2i (x) H 2 (x) + • • • + * n y (x) tf„ (x) 

V = 1 , 2 , •••,«), 

where •S'i,- (x) is the minor of <3? (x) with regard to the element in the ith. row 
and jth column, and H } - ( x ) is the function in § 3. Therefore, since [ <i> ( x ) | 
and all its minors are periodic functions, we obtain 

(-l)"+w » N / fc(x + l)H i (x+ 1) , \ 

^ (a;)= $(x) £*«<*H — H(*+D — +"v (j=l,2, "-' n) ' 

or by virtue of (15), 

*,(*) = (" D~|^«iOO + (- l)^|^ W2 (x) + ... 

*Py y ( X ) . <P n j (x) 

+ -*(^ (X) + •" + (_ 1} ¥(xT w " (a;) - 

Substitute the values just obtained f or 6\ ( x ) , • • • , 9 n ( x ) in (28) ; we find 
?>(z)=f^E(-l) m *i*(z)M*) 

-f^£(-i)* +1 *»(*)fo(*)+--- 

+ (- l)"" 1 !^ E (- l)^ 1 ^ (x) Z ifc (x) . 

If we make use now of the values of kj (x ) as given by (27), the last expression 
becomes, upon rearranging, 

- J/-.S Z^»y(a:)S(- l) m <Pift (z)*2* (a) 
+ 

+ (- D-^Tvf £M*) S (- D m w* (*)*.* (x) . 
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Since 

£(- l) w ^(*)»»(a) - |(^i)mi' (a;)) if i= z , 
the above expression for g%{x) easily simplifies^ and finally we have 

cji (x) = &>! (a;) A.-i («) + w 2 (a;) h i2 (*)+■•■ + co» (a; ) A in (a;) . 

Therefore j,- (a;) = jr,-i (a;) (i = 1 , 2, • • • , n) , so that we have the same 
solution obtained in § 3. 

This result can be stated as follows: When the equation (11) can be solved 
by means of a direct sum to the right, the solution of (1) obtained by the method of 
variation of constants is independent of the choice of the fundamental system of 
solutions of the associated homogeneous system, for in this case the series which 
occur are convergent, and therefore can be rearranged as above. 

§ 6. The asymptotic form of the solution in the left half plane. 

By an application of the results of the last section we can determine the 
asymptotic form of the first principal solution in the left half plane. 

Take as the solutions of (2) the intermediate solutions h'a ( x ) , h\j ( x ) , • • • , 
K j {x){j=\,2,---,n),oi%\. Put 

n 

(30) g'ii(x) = 'E,Z}j(x)h'i J (x) (*-l, 2, ••-, n), 

and apply the direct sum to the right in order to determine a[ ( x ) , • • • , 
w' n ( x ) ; a formal value of a'j (x) will be 

(3i) 5;( a) -(-iyH^ P<«+ 1) ^ ( * +1) +...l 

' L H'(x+l) J 

where H' (x) = | h'ij (x) | , and H'j (x) is the minor of H' (x) with respect 
to the element in the last row and jth column. The same convergence proof 
used in § 3 applies here, provided x is sufficiently distant from the real axis, 
since we have h' i} (x) ~ «# ( x ) with regard to v in the right half plane. 

According to the result of the last section this new solution must be identical 
with the first principal solution. As derived above it is known to exist only 
at a sufficient distance from the whole extent of the real axis. 

The typical term of (31) will again have the asymptotic form (16), the re- 
presentation holding with regard to x in the left, and with regard to v in the 
right half plane. To find the asymptotic form of the series, divide it into two 
parts: the first k terms, where k is such that x -{• k-\- 1 = Xiis a, point in the 
right half plane near the imaginary axis, and secondly, 4;he remaining terms. 

Take first the case where x approaches infinity along any ray making an 
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acute angle with the negative real axis. The first k terms can again be written 
in the form (20). By the definition of k we have in this case k < | x \ , so 
that the quantity 

! + ?) ( ™<*> 

which appears in the denominator of the typical term can be developed in the 
same way that the first term of (20) was treated, since all the expansions are 
valid if m<\ x\. Thus we have immediately that the first k terms of (31) 
when multiplied by h'ij (x) can be put in the form (21). 

Consider now the terms of (31) from the (&+ l)th on. They can be 
written 

y» i (xi-\- m) H'j (xi + m) 

m=o H' ( Xi + m ) 

Since Xi , Xi + 1 , • • • are points in the right half plane, the asymptotic form 
of the typical term of this series is again given by (16), the representation 
being now with regard to » . The series therefore takes the form 

« (e' L /p-) xl+m 

(32) 2~i ( I Wxt+mH-n-v M m \Xi-\- m ) , 

m—0 \ Xl I ' Il > J 

where M m ( Xi + m ) remains uniformly bounded for \v\ large. 

Let Xi = Mi + v V ' — 1, Tj = £+ fV— 1; then we can so choose Xi that 
jumi + £ — v > . We can also choose x so that the sign of \xv + f is the 
same as that of v . It can then easily be shown that the absolute value of the 
denominator of the (m+ 1 )th term of (32) is greater than \ Xl ^ +m ) +r >-" |. 
The absolute value of (32) is therefore less than 



M\(e*/ Pj y 

| x^ +r '- v | 



X? pj 



where M is a constant not depending on x . Since | e" / x f pj \ can be made less 
than unity, the sum of the last series is 

W\ m e I — 1 < M' p^+O^gx 111.] 

(33) 1-,.,1^+e- \\ Pi ) \ <Me \{pj \> 

where M' is bounded. 

The terms of (31) from the (k -{- l)th on are also to be multiplied 
by h'ij ( x ) , the absolute value of the principal part of which is 
I * l* u+f Kpje-^T \g-^^+0"". From (33) the absolute value of the product 
is seen to be less than M' | x |" u+f | pye - " | M_Ml . This vanishes more rapidly 
than any power of 1 / x along the rays considered, inasmuch as | pj e~ u \ is 
constant; for along these rays u increases negatively. Consequently the terms 
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of (31) from the ( k + 1 )th do not affect the asymptotic form of coj ( x ) h'a ( x ) , 
which is therefore given by (21).* 

It remains now to consider the case where x approaches infinity parallel to the 
negative real axis at a sufficient distance from it. A slight modification in the 
treatment of the first k terms must be made, for the quantity ( l+m/x) r ~ v+ ' ilx+m) 
cannot, in all the terms, be developed as before, for in the last terms | m / x \ 
will approach very close to unity. The modification to be made is apparent : 
we divide the first k terms into two parts; the first I of them, where, for instance, 
I is the smallest integer such that | x + I [ < ^ | x | (such a choice is always 
possible when x goes to infinity parallel to the negative real axis ), and secondly 
the remaining terms. 

The first I terms when multiplied by h'a (x) will again have the asymptotic 
form (21). We can show readily that the second group of terms when mul- 
tiplied by_the same quantity will contribute nothing to the asymptotic from 
of u'j (x) h'n ( x ) . The details of the argument are similar to the treatment 
of the terms after the kth, and we will merely outline them in a note.f 

For this case we therefore have again that Wy (x) h'jj (a;) is given by (21). 
From this it follows immediately that g'n(x) ~ g t (x) in the left half plane. 
Remembering that the new solution is identical with the first principal solution, 
we have: 

Theorem II. If /x > , the first principal solution is asymptotic to the formal 
series in the left half plane at a sufficient distance from the real axis.% 

*If f (x) = <p (x) + 6 (x) , where <p ( x ) ~ x" ( a + ai/x + • • • ) along some ray, 
while (x ) becomes small of order 1 / a'- x \ ( o > 1), then / ( x ) ~ x n ( a + ai / x + • • • ) 
since 1 / a' * I is smaller than any power of 1 / x for x large. 

t The terms in question can be written in the form (32), l^m^k , and x replacing 
xi . We can show first that the absolute value of the typical term would be increased if we 
changed (x + m)' ix + r — v in the denominator to (x + l)n x + r }-*, and the factor (x +m)> an 
to x$ m , where Xo is the smallest in absolute value of x +1, ••■ , x + k . We obtain then for a 
series greater term by term in absolute value, a geometrical progression, which can be summed, 
the sum being less than 

MIC"/*)--'! < MUxl2V-»*-rneH B Y\ 



'- — < 1 , and |x + /| < | x / 2 | . The product of this and the absolute value of the 
leading factor of the asymptotic form of h,j (x) is easily seen to be 

M|(x/2)^2^| = ^l f 

( x = u + v V — 1) . Since u is practically equal to — | x \ , and M is a fixed number, this 
decreases more rapidly than any power of 1 / x . 

t It should be noted that the asymptotic character of this solution of the non-homogeneous 
equation is somewhat different from that of the solutions of a homogeneous equation. For 
instance, the asymptotic form of T ( x ) holds only for — ir < arg'x < ir . That functions 
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It is interesting here to refer to (17). The quantities An (x) / A (x) 
( j , i = 1 , 2 , • • • , n ) , which occur there have zeros of order p at infinity. 
Consequently, if we assume that there exists a constant X such that lim x x gi{x) 
= ki, ki finite, as x approaches infinity parallel to the negative real axis, 
equation (17) gives 

iim^(^) = |:(-iy +i ^(i+^+---)iini^|^ 

+ lim(- l)"+ i+1 ^~ i (l + ^+ ■• ■)**+-*. 
; we must have 
X= n- v, k t = (- l)«+»+i- 



From this relation we see that we must have 

bA 



so that 

9 



a{x) = *-*[(- \Y +i+1 ^-+U{x)\ (*-l, 2, ••■, n), 



where f ; ( x ) tends to zero as | x | becomes infinite. From (24) and (25) it 
is clear this is obtained by breaking off the series gi(x) after the first term. 

§ 7. The second principal solution for p > . 

We shall now show the existence of a solution of (1) defined throughout 
the finite plane, and whose asymptotic form is known in the whole left half 
plane. It will of course differ from the first principal solution merely by a 
solution of the homogeneous system (2). 

Take as solutions of (2) the principal solutions on the left, that is, hy ( x ) , 
hij ( x ) > ' " • > h n j ( x ) ( j = 1 , 2 , • • • , n ) , where Ay ( x ) ~ »<,- ( x ) in the left 
half plane. Put 

n 

(34) g i2 (x) = ^Z}j{x)hij (x) (f = l,2, •••,»). 

We find in precisely the same way as before 

.b(x+l)H 3 (x+l) 



(35) Auj(x) = (- 1)"+' 



#(*+!) 



where H (x) is the determinant | Ay ( x ) | , while Hj (x) is the minor of this 
with respect to the element in the last row and jth column. 

If we attempt to solve (35) by the direct sum to the left, which we would 
naturally try in this case, the series obtained is found to diverge. We can, 

actually exist with the property specified in Theorem II can be shown by considering the 
series l/a; + l/(x + l) 2 + l/(a; + 2) 3 + '". It is easy to prove that this function is 
asymptotic to its formal expansion in powers of 1 / x , with regard to x , for — jr_< arg x < jt , 
and also for arg x = ± v , provided that | v I > o > 1 , where x = u + v V — 1 . 
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however, apply the integral formula (14). To prove this, we note first that 
the situation of the singularities of H } - ( x ) , and the fact that H ( x ) ~ | s# ( x ) | 
with regard to » in the right half plane, show that the expression on the right 
in (35) is analytic outside some region D, such as was described before. 
Accordingly the expression 

(36) *,(*) = (- !)•« C Ht+_l)H t (t+l) dt_ 

will formally satisfy (35), and in case the integral converges suitably it will 
represent an analytic solution. The contour L is the contour oo AB oo used 
earlier. 

To examine the nature of the integrand we make use of the intermediate 
solutions; for although, as noted above, H (t) is asymptotic to the deter- 
minant | an (t) | with regard to v in the right half plane, the minor Hj (t) 
is not asymptotic to the corresponding minor of | »« ( I > unless j = n, 
(see § 1). We therefore replace the elements of H, (x) by their values as 
given by (8). When we do this, the resulting determinant can be simplified so 
as to give finally 

(37) Hj (x) = H'j(x) + <p j+1 (x) H' J+l (x) + ■ ■ ■ + <p n (x) H' n (x) , 

where, as in the preceding section, Hj (x) is the minor of H ' ( x ) = | h'ij ( x ) 1 
with respect to the element in the last row and jth column, while <pj + i ( x ) , 
• • • , <p n (x) are periodic functions, and are analytic at a sufficient distance 
from the real axis. 

The first factor of the integrand in (36) can now be written 

(38) H!+±lM±±l) = ± Mt) Ht+_DHHt + i) 

When we make use of the asymptotic form of the quantities involved, the 
first term of (38) becomes, (see (16)), 



(39) 



(^/P,) m f J , 4 11 , r,j(t) 



(t+ l)*(* +1 >+ 



^(*+-J- + 5l t-). 



where v, (t) is limited, decreasing for large \v\(t—u+vl/— 1). When 
| v | is sufficiently great it is evident that (39) will decrease along Aoo and Boo 
in approximately the manner of an exponential function, inasmuch as we are 
assuming ju > . Since the functions <p i+ i ( t ) , • • • , <p n ( t) are periodic and 
limited along any line parallel to thereal axis, like results hold for all the terms 
of (38). The quantity ( e 2" (*-o ^-i _ i)-i re mains finite along L; conse- 
quently the integral will converge in such a manner that Wj ( x ) is an analytic 
function outside of the region D . 
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With these values of a } - (x) , equation (34) gives a solution of (1) analytic 
on the left outside of D . This solution can be extended to the right by means 
of the equations (1) themselves. From the fact that a# ( x ) and b ( x + 1 ) 
are analytic in the finite plane outside a certain circle, the extended solution 
can have singularities only at' points congruent on the right to the singularities 
of an (x) and b ( x + 1 ) . 

We now come to the asymptotic form of the solution. For this purpose we 
break up the contour L into two contours, L\ and L?. in the way that Birkhoff 
has done.* Here L\ is a fixed contour <x>AB°o, outside of which the integrand 
in (36) is analytic, and in case x lies between A<x> and B<x> extended to the 
left, L 2 is a loop circuit around x and the first point, x -4- I, to the left of AB , 
that is congruent to x; while if x is above Ax> or below jBoo , Z 2 is a loop circuit 
to infinity parallel to the positive real axis. In either case L\ and Z 2 are 
together equivalent to L , and L\ can be deformed so that no point congruent 
to x falls on it. 



00 

00 



x+l- 



L-, 



B m 

Fig. 3. 

Denote the parts of «/ ( x ) arising from L\ and Z 2 by Zfp (x) and w { P (x) 
respectively. Then 

n n _ 

(40) 9i2{x) = 'Ewf\x)h ij (x) + T l Z, ( p (x)hij(x) (i-l, 2, ••-,«). 



m=l H ( X + TO ) 



Consider first the case where x lies between A°o and 5oo; we have 

(41) 

The asymptotic form of the typical term is again given by (16). Also we see 
* The reason for this is evident. The singularities of the integrand in (36) are at the con- 
gruent points x , x + 1 , x + 2 , •••, due to the zeros of e 2T(*-'W-i — l ; and the fixed singu- 
larities of = '- , which are in the vicinity of .the positive real axis. The resi- 

H(t+1) 

dues at the first points will produce a series of the form (13), and the sum of the residues at 
the last points will be a periodic function. It is therefore natural to consider the two different 
groups of terms separately. 

Trans. Am. Math. Soc. 15 
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that precisely the treatment applied to the first k terms of (31) for the case 
where x approaches infinity parallel to the negative real axis, holds_for (41), 
provided AB is taken far enough to the left. Therefore w ( P (x)hi } - (a;) is 
again given by (21), so that we have 

n 

(42) £wf (x)li 3 (x) ~gi(x) <.i = 1,2, •••,«), 

3=1 

for x between A<*> and Boo. 

Next let x lie above Aoo or below Boo; then (41) is replaced by the infinite 
series 

apoo = ( - ipi f.iMi)tei). 

«»=i H( x + m ) 

Since the fundamental system Ay (x) , h^j (x) , • • • , h n j (x) can be expressed 
in terms of the solutions hy (x) , hy (x) , • • • , h nj (x) ( j = 1,2, •••,«), 
it follows from the results of § 5 that 

n _ n 

(43) Haf{x)h ij {x) = 1 ^03j{x)hij{x) = gi 1 (x) (i = 1, 2, •••, n). 

But we have found that gn (x) ~ gt (x) along any ray making an angle 
with the negative real axis, or parallel to that axis at a sufficient distance 
from it. We see then that (42) is true throughout the left half plane. 

We still have to examine the quantity up ( x ) . L\ is a fixed contour, so 
that Zfp (x) is periodic, and, moreover, its absolute value will have a finite 
upper limit. Furthermore, we can write 



-a), n = (-1)^' f b(t+l)H s (t+l) _dt 

°> KX) - ei-J-i J Zl H(t+1) e- 2 ^- 1 - e- 2 ^- 1 ' 



The integral approaches a constant, Cy , when a; approaches infinity along the 
negative imaginary axis, so that 

(44) lim e 2 ^~ l of} 3 ( x ) = C, (for arg z = ^) • 

From (36) itself we see 

(45) lim a ( j v (x) = Cj (for arg x = | J , 

where Cj is constant.* 

Consider now the product afp (x) hy (x) . It can easily _be shown that 
the leading factor £* x (pje~*) x x r > in the asymptotic form of ha (x) vanishes 
more rapidly than any power of 1 / x along any ray in the left half plane, includ- 
ing the positive axis of imaginaries. Then, since Z>P (x) remains bounded, 
uf-p (x) ha (x) will also vanish more rapidly than any power of 1 / x , and con- 

*The inversion of order of the double limits employed here is easily justified. 
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sequently it will not effect the asymptotic form of toy ( x ) ha ( x ) , so that 

ga(x) ~gi(x) (|^arg*<|Q. 

By using the second principal solutions of the associated homogeneous 
system, and expressing the sum formulae by contour integrals to the right, we 
have thus proved: 

Theorem III. If n > , there exists for (1) a second principal solution gu ( x ) , 
022 ( x ) , • • • , g n % ( x ) analytic throughout the finite plane, save for poles con- 
gruent on the right to the poles of a„ (x) and b (x -{- 1). This solution is 
such that ga (x) ~ gi (x) when x approaches infinity along any ray in the left 
half plane, inclusive of the positive imaginary axis. 

§ 8. The relation between the principal solutions. 

The difference between any two solutions of (1) is a solution of the homoge- 
neous system (2). To determine the relation between the principal solutions 
we make use again of the second principal solutions of the homogeneous 
system. We know that we must have 

n 

gt2(x) — gn(x) = ^8j(x) h tj (x) U = 1, 2, •••, «), 

where the functions B\ ( x ) , ■ • • , 8 n ( x ) are periodic. We have really already 
determined these periodic functions in the form of contour integrals; for a 
comparison of (40) and (43) with the above expression shows immediately that 
we must have 6 } - ( x ) = of 1 (x) when x is sufficiently distant from the real axis. 
The functions 0i ( x ) , • • • , 0« ( x ) are therefore analytic at a sufficient 
distance from the real axis, and as they have no essential singularities in a 
period strip, they are rational functions of e Ux -1 .* Thus we can write 

where my, n,-, nj, and qk are positive integers, and D$ , cPf are constants, 
while aikj is a fixed point. 

When we remember the identity of 6, ( x ) and w^ 1 ' ( x ) , for x sufficiently 
distant from the real axis, we see from (44) that we must have 

(46) mj 



'Osgood: Lehrbuch der Funktionentheorie, I, 1st edition, p. 404, or 2d edition, p. 466. 
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Likewi 


se from (45) we have 




df = 0, k<0, 


(47) 


X -1 X s kl 



[April 



0=1, 2, ••-, n). 



fei"(- c*'"*"'- 1 ) 1 C " 



The relation between the principal solutions is accordingly given by 
(48) 9ii (x) - ga (x) = ±1, (*)g t^^S^Ty 



with the conditions (46) and (47). 



(f = 1, 2, ••• , n), 



§ 9. The case p < . 

We have thus far explicitly assumed that ju is positive. We shall now in- 
dicate the treatment for /z negative. 

The series (15) diverges in this case, and so the direct sum to the right no 
longer furnishes a solution; but if we apply the direct sum to the left to 
(35), we can show, by using the asymptotic form of the terms, that the series 
so obtained will converge uniformly at a sufficient distance from the positive 
real axis. 

By means of a direct sum to the left we thus establish the existence of a 
solution analytic in the finite plane save for poles at the poles of a# ( x ) and 
b ( x + 1 ) and points congruent on the right. This solution corresponds to 
what was called the second principal solution when /z is positive. 

If now we take the first principal solutions as solutions of the homogeneous 
system, it can easily be shown that the resulting sum formulae can be taken 
to be contour integrals about a contour to the left. This gives a second 
solution of the non-homogeneous system analytic in the finite plane at a 
sufficient distance from the negative real axis, with its singularities at those 
of an ( x ) , b ( x + 1 ) , and the zeros of A ( x ) , or points congruent on the 
left. It thus corresponds to the first principal solution for ii > 0. 

These solutions will, as before, be asymptotic in the left and right half 
planes, respectively, to the series formally satisfying the equations, the first 
of them maintaining its asymptotic form in the right half plane at a sufficient 
distance from the positive real axis. The formal series have, however, a 
somewhat different form, and now are 
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gi(x) = Z y+ *(<7i + ^r+ •••) C*-l, 2, ...,n-l), 

For if we substitute expressions of this form in the equations, and compare 
coefficients of like powers of x, it is evident that at each stage we shall have 
n equations in which there are n unknown coefficients and coefficients pre- 
viously determined. Thus the coefficients can be calculated uniquely. In 
particular we have 

g%= ciinb ( i = 1 , 2 , • • • , n - 1 ), 

9=b. 

A case of interest, and one which will be used later in the treatment of the 
equation of the nth order, is that where an ( x ) = , for j =(= i + 1, i ! #= n. In 
that case the formal series become 

(49) gi(x) = i" fW Mf i + y+ •••• ) (* = i. 2, •••, n-l), 

where 

g% = fli, »+i fli+i, t+2 • • • flii-i, n & ( i = 1 , 2 , • • • , n), 

(50) 

0„ = 6 . 

We can also show the existence of solutions for the case ju < in a second 
way by referring it to the case n > . When we remember the definition of 
A (x) , Aji (x) , the system (17) can be written 

n 

Qi{x) = E««' (a0 Qj (x+ 1) + &i (x+ 1) (* = 1, 2, ••• , n), 
J=l 
where 

V'+i^ 
4 ' 



OLij (x) = X * ( a»y + ~ + ■ ■ ■ J , dij = ( — 1 ) 



-4»i& 



If now we put x = - (2+ 1), g { (x) = #* (- (z+ 1)) = h (z+ 1), 
and ju= — X(X>0), this system can be written 

71 

ff>(z+l) = T,c ij (z)g ] -(z) + b i (z+l) (i-l, 2, ..... n), 
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wher< 


e, when | z | > R , 








Ci y=(-iy +1 ^, 




b i ( ;z ) = z ^(b i + b ^+ ...), 


6,.= (-l)«+*i^ 



[April 



Since X is positive the previous existence theorems can be applied. The sub- 
stitution x = — ( 2 + 1 ) shows that the r61es of right and left half planes 
have been interchanged. 

We can summarize the above results as follows : 

When n < , first and second principal solutions exist and have their singu- 
larities at the same points as before. The first solution is asymptotic to the formal 
series in the right half plane; the second solution is asymptotic to the same series 
throughout the left half plane, and also in the right half plane at a sufficient distance 
from the real axis. 

■ § 10. The case n = . 

It still remains to consider the question of solutions for the intermediate 
case n = . It will appear that at times this case is exceptional. 

We shall first examine when the solutions, as found on the hypothesis that 
ix > 0, are valid when we put ju = 0. Consider the equation (15), and let 
us see under what conditions the series will still converge. The typical term 
can now by (16) be written 

(s+ m) y -> 

-^ [ dj-, x+ m ] . 

"j 

If | pj | > 1 , we shall have a series in which the terms decrease in practi- 
cally geometrical fashion, so that the series converges uniformly. Since 
|p» l< I P% l(* < «)■> the series for «i (x) , • • • , «„ (x) will all converge if 
| p„ | > 1 . The first principal solution as found for n > therefore still 
exists f or n = , if | p„ j > 1 . Under this same condition the second principal 
solution f or n > , given by the contour integrals, is also valid. For in this 
case (39) can be written 



■in /, ,dy Vj (t)\ 



and this expression decreases along A <x> and B oo in such a way that the 
integral in (36) converges. 
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We can then say: 7/ every root of the characteristic equation is greater than 
unity in absolute value, we can consider the case p = as a limiting case for p. >0 . 

In a similar way we can show that: When all the roots of the characteristic 
equation are less than unity in absolute value, the case p. = can be considered 
as a limiting case for p < . 

We still have to determine whether we can make the sum formulae valid 
if (1) some of the roots of the characteristic equation are greater in absolute 
value than unity and others less, or (2) some of the roots of the characteristic 
equation are equal to unity in absolute value. 

Consider the first case. Let pi , • • • , p& be greater than unity in absolute 
value, and pjt+i , • • • , p„ be less. Then, if we take as solutions of the homo- 
geneous equations the first principal solutions, we can use for the sum formulae 
for wi ( x ) , • • • , oik ( x ) the direct sum to the right. We will next show how to 
express o>k+i ( x ) , • • • , o> n ( x ) as contour integrals. As the work is similar 
for all of them, we shall consider the equation 

(5i) a^ (*) = (- D^ 6(a+ i ) ( fy 1 ( ; +1) . 

If x is in the right half plane and is sufficiently large, we know by (16) 
that the right-hand side of (51) can be written 

(a+ir*-' 

-=+1 [ dk+i; x + 1 J . 

Pk+l 

In § 1 it was stated that there is a set of intermediate solutions of (2) that are 
asymptotic to the formal series in the right half plane with regard to x, and 
with regard to v in the left half plane. Suppose these were the solutions of 
the homogeneous system that we were using; we could then write 

(51') Ae*n.i(*) = (- D^'^PWi *+ 1], 

where [dk+i; x + 1 ] is bounded for x large in the right half plane and for v 
large in the left half plane. Since | pk+i | < 1 , it is evident that for v large 
enough the above value of Ao>h+i (x) decreases, when x goes to infinity in a 
direction parallel to the negative real axis, approximately as an exponential. 
The equation could then be solved by means of a contour integral sum to 
the left. Return now to the consideration of equation (51). The first prin- 
cipal solutions can be expressed in terms of these intermediate solutions by 
equations similar to (8), and reference to (37) and (38) shows that the right 
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hand side of (51) would be replaced by a group of terms such as the one on 
the right in (51') in which every p is less in absolute value than unity. From 
this we see that cojt+i ( x ) can be expressed as a contour integral about a con- 
tour to the left, and similarly co^+2 (x) , ■ • ■ , w n ( x ) . 

If, therefore, we use the first principal solutions of the homogeneous system, 
and divide the resulting sum formula into two groups, according as they correspond 
to roots of the characteristic equation greater in absolute value than unity, or less, 
we can express the first group by a direct sum to the right, and the second by a 
contour integral to the left. Thus we have the first principal solution of the 
non-homogeneous system. 

In a similar way we can show that if we use the second principal solutions 
of the homogeneous system, and divide the sum formulae as before, we can 
solve the first group by a contour integral to the right, and the second by a direct 
sum to the left. This gives a second principal solution. 

The only remaining case is that in which some of the roots of the characteristic 
equation are equal to unity in absolute value. This case has not been com- 
pletely investigated for the homogeneous system. Here, in certain cases, at 
least, a modified sum may give a convergent series. An example of this is 
the equation occurring in the Weierstrass development of the gamma function, 
namely the equation 

g(x+ 1) - g (x) =-. 
The direct sum to the right is 

which diverges, while the modified sum 

i_i + (i__L\ + ... 

1 x^\2 x+lj^ 

is also a solution, and converges. 

When ju = the asymptotic form of the solution requires a slightly different 
treatment. Take the case where all the roots of the characteristic equation 
are greater than unity in absolute value, and consider the asymptotic form of 
the first principal solution. 

We need only find the asymptotic form of a 3 - (x) htj ( x ) , where 

^ (x + m)- r < f 

o)j(x) = 2- -^ [dj; x+ m\. 

m=l Pj 
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Take the first k terms of uj ( x ) , where % \ x\ — 1 < k <; \ \ x | ; we can write 
them in the form 

or 



,_r, k l / m \v-r 3 - 

-c- Z) -^ ( 1 + — ) [djix+m]. 

>,■ m=l Pj \ X ) 



In all these terms the expression (1 + m / x) v ~ rj can be developed in a series 
in 1/ x, since | m / x | < 1 . There is a marked difference, however, between 
these developments and those for the case p. > as given in § 4. In the 
present case each of the k terms of the summation when expanded will contain 
a constant term, while in the corresponding expression (20) in the former case, 
the first term alone when expanded contained a constant term, the others 
beginning with a; - ", • • • , ar (i_1)w . The final constant term in the present 

k 

instance will evidently be dj 2^ pj m • When k increases this approaches 

djl (pj — 1 ) , since \ pj \ > 1 . The first k terms of uj ( x ) , when multiplied 
by hij ( x ) , can therefore be written x" [ Xj* ; x ] , where Xj 4) approaches 
dj sn I (pj — 1 ) as A; increases. 

The terms of a, ( x ) after the kth behave ultimately like a geometrical pro- 
gression, since the ratio of successive terms obviously approaches 1 / pj as 
the number of the terms becomes large. It follows then that the remainder 
after k terms in the series for wy ( x ) , is, provided k is large enough, of essenti- 
ally the order of the ( k + 1 )th term, which is M ( x + k + 1 ) y - r VpJ +4+1 , where 
M is bounded. The product of this and hj(x) can be written 



,-, *" (. . k+\\- r > 



where | M \ is less than some fixed constant. Since k increases in a constant 
ratio to x , this becomes smaller than any power of 1 / x , so that we have ulti- 
mately coj(x) hi, (x) = x" [\j h) ; x] where Xf approaches djSij/ (py — 1) for 
k large. 

In this way we can determine the asymptotic form of the solution, and we 
are led to try to find formal solutions of the form 






(i = 1, 2, ••• , n). 



That the coefficients in these expressions can be uniquely determined so as 

to satisfy the equations can be shown easily by substitution. In particular 

we find 

bA- 
(52) g t = (- l)»+ i+1 -^ (i-l, 2, .-., n), 

A 
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where A is the determinant | a# — 5# | , ( 5,-,- = , if i #= j , 8« = 1 ) , and Ai 
is the minor of A with reference to the element in the last row and ith. column. 
Also we find again that the first and second principal solutions are asymptotic, 
in the right and left half planes respectively, to the series formally satisfying the 
equations. 

§ 11. The equation of the nth order. 

We shall consider a single non-homogeneous equation of the nth order, 
which we write in the form 

(53) g (x + n) + a! (x) g (x + n — 1)+ ••• + a n (x) g (x) = b (x) . 

The functions ai (x) , • • • , a„ (x) , b (x) are assumed to be rational, so that 
we have 

■ / . e (I) 
at 



( a (1) \ 

i i (x) = x i »^a i + ^- + •••J, 

b(x) = x" ^6 + —+•••), 



(for \x\ > R) . 



In place of the single equation (53), consider the system 
gi (x+ 1) = x» g 2 (x), 
2 (x+ 1) = x* g 3 (x), 



. 1 (x+l) = z*gA*), 



z»(a n _ m + %+•••) 
(54) fl-„(a;+l) Ey-; • jw ^r 7 nujfW 

{( i+ -T i )( i+ i)-( i+ V)} 

/ 6 a > \ 

j.v-c-W 6 + —4 ) 



If from these equations we eliminate g 2 (x) , • • • , g n (x) , we shall find that 
gi (x) satisfies (53), so that the system is equivalent to the original equation of 
the nth order. When we expand the quantities ( 1 +j/x )" , which occur in 
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the last equation above we find that that equation can be written 

g n (x+ 1) = ai(x)gi(x) + d 2 (x)g 2 (x) + h a n {x)g n (x) + b (x + 1) , 

where 

di(x) = — x» ( a„_ i+ i + ~ + • • • J > 



( 6 (1) \ 



so that the system is of precisely the form (1). The determinantal charac- 
teristic equation reduces to 

(55) p n + a lP n ~ 1 + ■•• + a n _! p + a* = 0, 

which can therefore be found directly from (53). 

We can now apply to the system the existence theorems already estab- 
lished, and we thus have solutions g\ ( x ) of (53) in all cases except when n is 
zero and some of the roots of the characteristic equation are equal to unity in 
absolute value. 

The asymptotic form of the solutions of (54) can be determined by applying 
the general results for a system (for n < 0, apply (49)). In each case it is 
found that the solution gi(x) is asymptotic to the series formally satisfy- 
ing the equation (53). This formal solution is now 

g(x) = x v - n »( g + — + ••• J (for ,«>0), 

g(x) = x v (g + 9 —+ •••) (for„5i0) ( 



where 



b 
g = — ( f or M > ) , 



a n 



g = b ( f or n < ) 

h 

( for /x = ) . 



l + ai + a-i+ ■■■ + a n 



The constant term can in all cases be found either from the known form of 
the constant term of the solution of the corresponding system (for n < 0, 
apply (50); for n = 0, use (52)), or by direct substitution in equation (53). 
We can now state: 
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Theorem IV. If fx < , and all the roots of the characteristic equation (55) are 
different from zero, the equation (53) has two principal solutions gi(x) and g[ ( x ) . 
The solution gi (x) is analytic in the finite plane save for possible poles at the 
poles of a\ ( x ) , • ■ • , a n -\ ( x ) , b ( x ) , at the zeros of a„ ( x ) , or at points con- 
gruent on the left; while g[ (x) is analytic save for poles at the poles of ai (x) , 
• • • , a n (x) , b (x) or points congruent on the right. If n = , the same is 
true provided in addition that none of the roots of the characteristic equation are 
equal to unity in absolute value. There exists a series g ( x ) formally satisfying 
(53), and in all cases in which we have shown the existence of solutions we have 
gi (x) ~ g (x) in the right half plane, while g\ ( x ) ~ g ( x ) in the left half 
plane. 

Princeton University, 
May, 1912. 



